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We show that several aspects of the low-temperature hydrodynamics of a discrete Gross-Pitaevskii 
equation (GPE) can be understood by mapping it to a nonlinear version of fluctuating hydrodynam¬ 
ics. This is achieved by first writing the GPE in a hydrodynamic form of a continuity and an Euler 
equation. Respecting conservation laws, dissipation and noise due to the system’s chaos are added, 
thus giving us a nonlinear stochastic field theory in general and the Kardar-Parisi-Zhang (KPZ) 
equation in our particular case. This mapping to KPZ is benchmarked against exact Hamiltonian 
numerics on discrete GPE by investigating the non-zero temperature dynamical structure factor and 
its scaling form and exponent. Given the ubiquity of the Gross-Pitaevskii equation (a.k.a. nonlin¬ 
ear Schrodinger equation), ranging from nonlinear optics to cold gases, we expect this remarkable 
mapping to the KPZ equation to be of paramount importance and far reaching consequences. 


Introduction: Low dimensional classical and quantum 
systems are often very counter-intuitive and different 
from their higher dimensional counterparts [T] . One such 
example is the width of the line-shape of the phonon 
peaks in the dynamical structure factor. Contrary to 
the expected k‘^ behavior in higher dimensions [2], the 
power is anomalous in low dimensions. Linearized hy¬ 
drodynamics, which predicts a diffusive broadening, fails 
in one dimension (ID). This immediately creates a need 
for a nonlinear hydrodynamics that could describe low¬ 
dimensional systems. Such a theory beyond the conven¬ 
tional Luttinger Liquid would describe the super-diffusive 
broadening in low dimensional systems. 

An experimentally realized [3] system that provides a 
remarkable platform for probing low-dimensional fluids 
is the system of a ID weakly interacting Bose gas at 
non-zero temperature. Using a variant of Bragg spec¬ 
troscopy SIS] one could probe the dynamical structure 
factor of the Bose gas, thereby unraveling the nonlin¬ 
ear phenomenon in low-dimensional fluids. Needless to 
say, the underlying theory that describes [6] this cold 
atomic system, namely, the Gross-Pitaevskii equation 
(GPE) or the nonlinear Schrodinger (NLS) equation is 
ubiquitous in areas such as optics, cold gases and math¬ 
ematical physics. Although the strictly continuum GPE 
is integrable, the experimental realizations break integra- 
bility in one or more ways such as, presence of a lattice or 
trapping potential, energy loss, escape or unwanted evap¬ 
oration of particles. Here we focus in particular on the 
discrete (lattice) version of GPE which is not integrable; 
such a discrete GPE has been realized in experiments on 
waveguide lattices [?]• 

The ubiquity of such equations and cutting edge 
technologies available to probe statistical properties of 
such systems, enhances an urgent need for writing a 


stochastic nonlinear theory that makes transparent the 
role of various components that result in a complex 
nonlinear-driven-dissipative phenomenology. Establish¬ 
ing this strong connection between GPE and stochastic 
nonlinear differential equations (which turns out to be 
a two-component KPZ equation in our case) helps in 
using the tools available in the literature to make far 
reaching predictions about the statistical mechanics of 
systems such as a ID Bose gas or optical waveguides. 
In the converse, one could also use such systems as an 
experimental test bed for KPZ phenomena, providing 
much needed additional experimental realizations of 
KPZ physics [HHin]. 

In this Letter, we analyze the low-temperature hy¬ 
drodynamics of the GPE, which is known to be a 
valid description for systems such as a ID weakly 
interacting Bose gas or optical waveguides. We present 
a discrete GPE that governs the dynamics of such 
complex wave fields (which are atomic fields in the case 
of cold atoms or optical fields in waveguides). We write 
down continuity-like and Euler-like equation for the 
macroscopic density and velocity fields and derive the 
nonlinear fluctuating hydrodynamics. The coefficients of 
the resulting nonlinear fluctuating theory are expressed 
in terms of the underlying parameters of the system 
(such as coupling strength, background density). Having 
established this, we present results for the dynamical 
structure factor S{k,uj) (i.e., fourier transform of cor¬ 
relation function of fields obeying nonlinear fluctuating 
hydrodynamic theory), namely its scaling function and 
the underlying anomalous exponent. This effective 
nonlinear hydrodynamic theory is finally benchmarked 
against exact Hamiltonian numerics, which also supports 
a recent remarkable conjecture that the long-wavelength 
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dynamics of a classical ID fluid at finite temperature is 
in the Kardar-Parisi-Zhang (KPZ) universality class [TT]. 
In addition to the confirmation of the 3/2 exponent, we 
have taken a big step forward in showing agreement with 
the Prahofer-Spohn scaling function [12]. Therefore, the 
notoriously difficult problem of computing the dynamical 
structure factor (or density-density correlations) can 
now be connected to correlation functions of familiar 
stochastic differential equations. We discuss certain 
points such as role of integrability, analysis of Lyapunov 
exponents, and future challenges and present the conse¬ 
quences of this mapping and its possible exploitation in 
understanding experiments governed by the Nonlinear 
Schrodinger equation. 

Nonlinear fluctuating hydrodynamics and GPE: The 
semi-classical Hamiltonian describing a strictly one¬ 
dimensional gas of bosons of mass m and contact 
interaction strength g is given by 


H = 



Idx'tpf 

2m 




(1) 


which in conjugation with Poisson brackets 
{'0*(x),'0(^)1 = i5{x — y) gives the time-dependent 
GPE, 


idt^ 


1 

2m 


dl'>P■ 


( 2 ) 


This is an integrable system. But all physical realiza¬ 
tions are not in this ideal limit: they may be in a lattice 
rather than with continuous translational invariance, and 
they are not strictly one-dimensional, so the interaction 
has a nonzero range. Here we will assume we are not 
in the ideal integrable limit, so integrability is destroyed 
and this nonlinear classical system is chaotic at nonzero 
temperature. The specific integrability breaking we con¬ 
sider is the discrete GPE (equivalently, NLS) on a one¬ 
dimensional lattice (see Eq. but the results should 
apply more generally. Eor optical applications, ^ is a Kerr 
nonlinearity and is the intensity of the light field. 

We examine the hydrodynamics of the equilibrium 
steady-state that this chaotic system approaches at long 
times. We are interested in the hydrodynamic scaling of 
the density-density correlation 


s{x,t) = mx,t)f\mo)\^) - {\mo)\y ( 3 ) 


with (•) denoting the average over the statistical steady 
state. '0(x, t) = p(x, defines the density p(x, t) 

and the phase 0. The velocity is v{xN) = ^ • We 

work at low enough temperature that the rate at which 
phase slips occur at equilibrium is negligible, so the ve¬ 
locity is a conserved quantity, as is the density. The 
continuity and Euler equations are 

dtp + dx{pv) = 0, dfV + = 0. (4) 


The equilibrium state has average density po = (|'0(x)p) 
and we consider the case of zero average velocity. 

In the regime we are considering, Eq. (© refers to 
small deviations from the average density. Hence if we 
linearize Eq. Q, taking p ^ po g and ^ 0 + i;, we 
obtain 


dfU + dxAu 


with u = 




Eq. (§ gives right- and left-moving sound modes with 
speed c = \fgpoJm. One can view Eq. ^ as the dynam¬ 
ics of a linearized Luttinger liquid whose S{k,uj) consists 
only of a pair of delta function peaks at cj = ±c|/c|, cor¬ 
responding to undamped phonons. We need to add to 
this, the scattering between the phonons due to the non- 
linearities. In linear fluctuating hydrodynamics one adds 
damping and noise to Eq. ^ which broadens the sound 
peaks in S{k,uj)^ giving them a line width that scales 
“diffusively” as T{k) ^ k^. This works fine in dimension 
d > 3 [13], but fails in one dimension [14]. An example 
showing this anomaly is simulations of Eermi-Pasta-Ulam 
(EPU) chains which report super diffusive broadening of 
the sound peaks [ISHIS]- A EPU chain consists of masses 
coupled to their nearest neighbors through anharmonic 
potentials. The discrete GPE has a similar structure, 
although instead of being anharmonic in the displace¬ 
ments, it is anharmonic in the local density. To capture 
such anharmonic behavior, it has been proposed recently 
to use a nonlinear extension of fluctuating hydrodynam¬ 
ics [20]. We will follow this strategy to obtain the hy¬ 
drodynamic scaling of S{k,uj), which we then compare 
to exact Hamiltonian numerics. The prescription of non¬ 
linear fluctuating hydrodynamics [20] consists of adding 
diffusion and noise matrices in Eq. [^giving. 


dfU + dx 


Au + ^ ^ dx i^Du^ -f- 


CX,P = 1 


= 0 


(6) 

where D and B are diffusion and noise matrices. Above, 
the vector = du^du^p with j = {gv, 

Dropping the quadratic terms would correspond to lin¬ 
ear fluctuating hydrodynamics and thus yield diffusive 
sound peaks. Eor our application we are interested in the 
stationary, mean zero process governed by Eq. (H , again 
denoted by (p(x, t),'r’(x, t)). The equal time, static cor¬ 
relations are expected to have short range correlations. 
Hence we define the susceptibilities 


Cl = J dx{{g{x,0)g{0,0))-pQ) ,C2 = J dx{v{x,0)v{0,0)) 

( 7 ) 

where the cross terms vanish because g{xN) and v{xN) 
have different parity. The fluctuation dissipation relation 
is given by DC CD = BB^ {C is a diagonal matrix 
containing Eq. [^. In addition, space-time stationarity 
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enforces in general the relation AC = which im- 

2 

plies C 2 = A Cl. In Eq. (4) the linear terms dominate and 

Po LI 

to obtain better insight to the solution one has to trans¬ 
form to normal modes which have a definite propagation 
velocity. Thus one introduces a linear transformation in 
component space, by setting 


(«) 

such that R satisfies RAR~^ = diag(—c, c). In addition 

we require that the ^-susceptibilities are normalized to 

one, which means i?diag(ci, ^ci)RA = 1. Up to an over¬ 
do 

all sign, R is uniquely determined and given by 


R = 


1 

c^/2^ 




( 9 ) 


Then the equation for the normal modes (i.e, the left and 
right chiral sectors) reads 


dfC^cr + dx [crC(/)cr + — dx{Drot<P)a + (^rotOcr] — 0 

( 10 ) 

with cr = ± and indicating the matrices rotated by 
R matrix, Dr-ot = RDR and 5rot = FIB. The coupling 
matrix is given by. 




2po 


3 1 

1 -1 


r + = -— / L 


-1 1 
1 3 


( 11 ) 

Since Eq. (10) is nonlinear, it is still difficult to com¬ 
pute the covariance ( 0 cr (^5 O0cr'(O, 0)) for the mean zero, 
stationary process. Although obvious, one central obser¬ 
vation is that in leading order the two peaks in S{xA) 
separate linearly in time. Hence, in the equation ( pT] ) 
for (/)cr, the terms and (0-cr)^ tarn out to be ir¬ 

relevant [19] compared to (0cr)^- Albeit they may effec¬ 
tively renormalize the non-universal coefficients (in front 
of all terms), they do not impact the universal proper¬ 
ties. Therefore, preserving universality we can decouple 
Eq. [T^into two components giving. 


dt<pa + dx [cTC^cr + ~ (U^rot0)cr + (^rotOcr] “ ^ 

( 12 ) 

Eq. is the stochastic Burgers equation (KPZ in 
“height function” where ha- = Oxc/^a) and for it the 
exact scaling function is available. 


{(paix,t)(pa{0,0)) = i>^t) ^/^/KPz[(Ai) ‘^^^{x-act)] 

(13) 

valid for large xR. A is a non-universal coefficient, which 
here is explicitly calculated to be 


A = 2^/2|G^J (14) 


One should keep in mind that the value of A derived 
above (0 will get renormalized m due to the dis¬ 
carded non-linearities as explained above. Note that A 


does not depend on I^rot or ^rot- This says that, while 
some dissipation and noise is needed to maintain sta- 
tionarity, the asymptotic form of the correlation is dom¬ 
inated by The definition of /kpz is somewhat in¬ 

direct. One first computes a family of probability densi¬ 
ties Px{s) indexed by x, in essence defined by a Ered- 
holm determinant which has to be evaluated numeri¬ 
cally [29]. Then /kpz(^) = f dspx{s)s‘^ . /kpz has 
the following properties: /kpz ^ 0, / dxf kpz{x) = 1, 
/KPz(a:) = fKPz{-x), f dxfKPz{x)x'^ = 0.510523 and so 
on. /kpz looks like a Gaussian, but with a large \x\ 
decay instead as exp[—0.295|xp] [12]. Eor a few dis¬ 
crete models there is in fact a proof [H mi- Eor the 
stochastic Burgers equation there is a tricky replica com¬ 
putation yielding the same result [22]. In some molecu¬ 
lar dynamics simulations, one records directly the cor¬ 
relation S{xA) [23] where x is the space coordinate. 
But more conventionally, as also relavant in this Let¬ 
ter, one studies the structure function, which is defined 
as the space-time Eourier transform of S{xA). We de¬ 
fine S{kA) = S{xA)dx. As argued before, the 

asymptotic scaling form is expected to be of the form 

S{k,t) = A(e-'''=‘+e-‘'''=*)ci/KPz(fc(A|t|)2/3). (15) 

The two sound peaks are symmetric reflections of each 
other. Considering only the right mover and setting = 
ck, one concludes 

S{k,uj +UJ],) = J dte“^qci/KPz(fc(A|ty/^) (16) 

= I dte‘(“/"l'=l^'^)*(A|fc|3/2)-llci/KPz(|t|^/3). 

Thus defining 

h{uj) = y dte'“‘/KPz(|ip^^), (17) 

one arrives at 


S{k,uj + cjk) = |ci(A|/c|^/^) ^h{uj/X\k \^^‘^). (18) 


If the maximum of S is normalized to 1, then the prefac¬ 


tor in (18) is set to 1 and h is replaced by h/h(0). 

One should cautiously note that the decoupling hy¬ 
pothesis is a subtle issue. The (j)a{xA) fields fluctuate 
without any spatial decay. It is in fact, only the correla¬ 
tions that are peaked near ±ct. However the decoupling 
of the components can be seen directly on the level of 
mode-coupling in the one-loop approximation. As sup¬ 
ported by numerical solutions [23], it is safe to use the 
diagonal approximation 


{(j>a{x,t)(pa'{^,Q)) = daa' fa{x,t) . (19) 

Of course fa{x,t) = f-a{—x,t). In one-loop one has {p 
is the phenomenologically added dissipation) 

dtfa{x,t) = {-acdx + \vdl)fa{xR) 

I J dyf„{x -y,s)dlM{y,s) (20) 


+ 
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FIG. 1: (Top) Comparision between exact numerics of dis¬ 
crete GPE with nonlinear fluctuating hydrodynamics. Here 
the temperature T = 0.005, L = 5 x 2^^. The best fit to the 
scaling function h{(ju)/h{0) is given by Xopt ~ 0.005 (and theo¬ 
retical A ~ 0.045) with universal tail (shown by Log-Log inset) 
close to Ct;“3 at significantly large frequencies p^ . (Bottom) 
Similar comparison for a different set of parameters. Here 
the temperature T = 0.001, L = 10 x 2^^. The best fit to the 
scaling function h{(ju)/h{0) is given by Xopt ~ 0.0028 (and the¬ 
oretical A ~ 0.0041). The wave vectors are given by /c = 
where the values of the integers q are given in the legends. 


with the memory kernel 

M{x,t)=2 ^ (21) 


crcr' = ± 


The terms with a ^ a' have a very small overlap. But 
the diagonal terms proportional to (GJ^r)^ contribute 
to the long time behavior. By explicit computation one 
checks that the self-interaction term dominates the mu¬ 


tual one. Eq. (20) can be studied numerically by an iter¬ 


ation scheme. The asymptotic shape of the sound peak 


is, of course, not the true scaling function /kpz, but so 
with a relative error of about 4% [23] . It is of utmost im¬ 


portance to have such a deterministic expression (20) for 


the correlators of ID Bose gas that captures the physics 
beyond a conventional Luttinger Liquid. This computa¬ 
tionally advantageous method (Eq. [2Q| ) along with our 
derived Eq. [T^and the above established stochastic non¬ 
linear field theoretic description of ID Bose gas needs to 
be benchmarked against brute-force Hamiltonian numer¬ 
ics of the underlying GPE Hamiltonian. 

Hamiltonian numerics of discrete-GPE: We now go to 
the discrete version of above time-dependent GPE (Eq. 

that now governs the dynamics of a complex-valued 
'0(n,t), with integer n = and periodic bound¬ 

ary conditions. Discretization is achieved by substituting 
X ^ na where a is the lattice spacing and Na is the sys¬ 
tem size L. The discrete version of time-dependent GPE 
reads, 






+ g\'ip{n,t)\‘^ij{n,t) (22) 


where kq = j^q for integer q = ( — y + l),^ and 
XE denotes our inverse-Fourier transform, XE{G{kq)} = 
b ^ (slightly unconventional due to 

explicit presence of lattice spacing a). The local energy 
and the local number density, |'0(n)p, are conserved. Ac¬ 


cording to standard classifications, Eq. (22) is listed as 
not integrable [24|. Hence one would expect that H 
and N = \ are the only conserved fields 

and that the set of equilibrium states is of the form 
^ > 0, /i G M, in the limit of large N. 
Therefore, the above discretization scheme for the inte¬ 
grable continuum GPE breaks the underlying integrabil- 
ity. In fact, we find that, in order to make connection 
to fluctuating hydrodynamics and subsequently KPZ, we 
require broken integrability and the resulting chaos. 

In this section, we describe the Hamiltonian exact nu¬ 
merics m starting from Eq. 


22 The time evolution 


is obtained by the well-known leap frog splitting tech¬ 
nique where the system is evolved alternatively (setting 
^ = m = 1, and always choosing po = 1) by kinetic, 

and potential, 

t), terms in sequence Vx • 7^ • Vx, with 

time step r. 

In the simulation [25] we measure the structure func¬ 
tion S{kq^uj): At each time step we obtain the time 
evolved density p{nH) = which we then space- 

time Fourier transform to p{kq^uj). Then the dynami¬ 
cal structure factor is the ensemble average: S{kq^uj) = 
{\p{kq,uj)\‘^). These results are expected to depend only 
on the total energy and particle number of the initial 
condition, since the dynamics are chaotic and we expect 
ergodic. The chaos for our parameters a = 5 and r = 2 
has been confirmed by observing positive Lyapunov ex¬ 
ponents. 
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For the random initial conditions we assume that the 
Fourier coefficients Qk^Ok are independent Gaussian ran¬ 
dom variables with mean 0 and covariance given by 






1 


2pqL Cik 


(23) 


^+P0 


where = \j\ {^ + 2/9o) and ak^ = ^ 

In Ref. [25] the low temperature dynamical structure 
factor was simulated numerically and the KPZ scaling 
exponent was observed, i.e. phonon line width Tk ^ \k\^ 
with ^ = 1.510 ±0.018. Here we provide a more quan¬ 
titative comparison with the full scaling function (Fig. 
[^. We also have fully outlined the mapping to nonlin¬ 
ear fluctuating hydrodynamics, which tells us that the 


structure factor should be of the form shown in Eq. (18). 


This means that the structure factor we obtain from our 
simulations must scale with the KPZ scaling exponent 
and scaling function in the hydrodynamic limit. 

In Fig. I^we show the remarkable quantitative agree¬ 
ment between exact Hamiltonian numerics and the ex¬ 
pectations of a nonlinear hydrodynamic theory with fluc¬ 
tuations. Our results are in a regime where the system 
is not near integrability, due to the large lattice spac¬ 
ing a. We have checked that on approaching integrabil¬ 
ity by reducing a we find strong deviations from KPZ, 
both in terms of scaling form and exponent, as we enter 
the regime of crossover between integrability and KPZ 
scaling. The discrepancy between the optimally chosen 
value of A (Xopt) and the one expected from KPZ cor¬ 
respondence (Eq. 14) probably arises due the fact that 
higher-order nonlinearities and the different chiral sectors 
effectively renormalize the first relevant nonlinearity of 
the specific chiral sector under consideration. Such a dis¬ 
agreement has also been seen recently in case of the EPU 
problem m- One does require an effective renormaliza¬ 
tion scheme to make more precise connections between A 
and Xopt • 

In conclusion, we demonstrate a strong connection 
between the statistical mechanics of a discrete NLS/GPE 
and a nonlinear hydrodynamic theory with fluctuations. 
This was done by first formulating the GPE in terms 
of hydrodynamic variables (conjugate classical fields) 
and then adapting a recent procedure in formulating 
a fluctuating version of the nonlinear hydrodynamic 
theory |20]. In our case, the resulting theory is shown 
to be of the KPZ universality class. This immediately 
enables us to use the rich physics of KPZ class and a 
well-established one-loop approximation to make pre¬ 
dictions for GPE. This was then benchmarked by exact 
numerics of the molecular dynamics type. Given the 
wide range of phenomena described by these equations, 
our results have implications in fields ranging from cold 
gases to nonlinear optics. Moreover, extending this 
mapping to coupled Nonlinear Schrodinger equations 


(also an experimentally realized situation in both cold 
atoms na HI] and nonlinear optics [28]) is shown to 
give a whole zoo of interesting dynamical critical phe¬ 
nomenon arising due to coupled stochastic differential 
equations [3T] . 
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